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^ ■ Abstract 



We present an analysis of the dependence of the resonant-orbital order and magnetic scattering 



'-^ ' spectra on the spin configuration. We consider an arbitrary spin direction with respect to the local 

c 

O ' crystal field axis, thus lowering significantly the local symmetry. To evaluate the atomic scattering 

o : 

in this case, we generalized the Hannon- Trammel formula and implemented it inside the framework 

'. 

, of atomic multiplet calculations in a crystal field. For an illustration, we calculate the magnetic and 

G\ 

Q^ ' orbital scattering in the CE phase of Lao sSri 5M11O4 in the cases when the spins are aligned with 

\o ' 

■ the crystal lattice vector a (or equivalently b) and when they are rotated in the afr-plane by 45° with 

O ; 

. respect to this axis. Magnetic spectra differ for the two cases. For the orbital scattering, we show 

-i— > 



that for the former configuration there is a non negligible 0^0' (n — > 7r') scattering component, 
which vanishes in the 45° case, while the a — ► ir' (ir — > a') components are similar in the two 
cases. From the consideration of two 90° spin canted structures, we conclude there is a significant 



q ', dependence of the orbital scattering spectra on the spin arrangement. Recent experiments detected 



a sudden decrease of the orbital scattering intensity upon increasing the temperature above the 
■ Neel temperature in Lao.sSri.sMnO^ We discuss this behavior considering the effect of different 

types of misorientations of the spins on the orbital scattering spectrum. 

PACS numbers: 61.10.-i,75.10.-b,75.25.+z,71.27.+a 



I. INTRODUCTION 



The physics of perovskite-type manganites is one of the very exciting areas of the present- 
day solid-state research, displaying abundance of interesting phenomena, of which the most 
known is the huge resistivity change upon the application of magnetic fields, the so-called 
colossal magnetoresistance effect^. Manganites have rich phase diagrams with several types 
of magnetic, charge and orbitally ordered phases, often exhibiting mesoscale phase coexis- 
tence. The physics of these systems is controlled by a complex interplay of magnetic, orbital, 
charge and structural degrees of freedom. Orbital order is characterized by the spatial ar- 
rangement of the outermost valence electrons. Although it was predicted and described 
theoretically 50 years ago^, it has been recognized as a hidden degree of freedom, because 
of difficulties related to its direct observation. 

The resonant x-ray scattering (RXS) is a promising tool for the study of orbital order, and, 
at the same time, convenient for the study of magnetic and charge ordering. In particular, 
it can be applied to investigate the correlations between the orbital order parameter and 
the magnetic and charge order parameters. The first RXS experiments with the goal of 
detecting orbital ordering were done at the K-edge. At the manganese K-edge, one Is 
electron is excited into the empty Ap band after absorption of a photon of the resonant 
energy. This intermediate state has a finite lifetime, after which the ion emits a photon and 
relaxes back to the initial state. Today it is believed that the observed signal is not sensitive 
to the onsite orbital ordering of the 3c? orbitals, but to band-structure effects including 
the effect of the Jahn- Teller crystal field on the 4p states together with the hybridization 
of the 4p orbitals with the neighboring oxygen orbitals (also controlled by the Jahn- Teller 
distortion) and with empty 3d orbitals of the neighboring Mn atom o?i 4 ^ . Instead, the L/2,3- 
edge scattering measurements ^ i^Pi 1 ^ 11 , stimulated by the theoretical work of Castleton and 
Altarelli^, are sensitive to both the onsite orbital ordering and the Jahn- Teller distortion. 
At the Mn L 2j 3 edges, the virtual dipole transition consists of an excitation of one 2p electron 
into the 3d shell, thus directly probing the orbital ordering. The specific contributions of the 
orbital ordering and Jahn- Teller distortions must then be disentangled in order to interpret 
the results of such experiments. This is usually done by comparing the experimental spectra 
to theoretical simulated spectra. Previous theoretical wor k 1 ^ 11 has indicated that the Jahn- 
Teller distortion affects almost exclusively the L 3 edge. However, to monitor the effect of 



a change in the magnetic order on the orbital ordering in particular (or on the Jahn- Teller 
distortion), one needs also to know the direct influence of the magnetic configuration on the 
orbital spectra. So far, to our knowledge, this effect on the L edge orbital scattering has 
not been addressed. Also, obviously, if the spin configuration has a significant influence on 
the orbital scattering, it must be properly included in the analysis used to disentangle the 
effects of the orbital ordering and Jahn- Teller distortion. 

Recently, Wilkins et a/.— reported measurements of the temperature dependence of the 
orbital scattering at the Mn 1^,3 edge in Lao.sSri.sMnC^ which displayed a substantial 
decrease in the orbital scattering intensity, as the temperature is raised, at the transition 
from the antiferromagnetic to paramagnetic phase. Analogous behavior when crossing the 
Neel temperature was previously reported in LaMnO^, at the Mn K edge, and in KCuF s 14 ^ 
at the Cu K edge (but the origin may not be the same). More recently, Staub et a/.— also 
presented measurements of the temperature dependence of the orbital scattering at the 
Mn L/2,3 edge in Lao.sSri.sMnC^, which showed a smaller change at the Neel temperature. 
Furthermore, based on a polarization analysis, they showed that the orbital signal below the 
Neel temperature includes a non negligible contribution of magnetic scattering from twinned 
minority ferromagnetic domains. However, it is still unclear whether this fully explains the 
temperature dependence of the observed orbital signal*^. In principle, the change in the 
magnetic configuration at the Neel temperature in the majority domain could also affect the 
spectrum. 

In this work, we address the direct influence of the magnetic configuration on the orbital 
scattering spectra for a fixed orbitally ordered configuration. As a prototype system, we focus 
on Lao.5Sri.5Mn04. We find that the spin configuration has indeed an influence on the orbital 
scattering, but mostly when the spins are noncollinear. For the magnetic scattering we find 
that the spin configuration has a non-negligible influence on the spectrum, and this even 
in the collinear CE structure. Previous studies generally considered spin oriented along the 
highest local symmetry axis of the Mn 3+ crystallographic site. In order to allow an arbitrary 
orientation of the spin with respect to the local crystal field axis, we generalized the Hannon- 
Trammell formula^ for atomic scattering, derived for spherical local site symmetry, by a 
systematic expansion in invariants of decreasing symmetry down to Q. We use these results 
to derive the scattering formulae for different spin configurations, including the polarization 
dependence. Finally, we consider various types of spin misalignments and misorientations to 



understand the effect of spin disorder on the orbital spectrum of the paramagnetic phase. We 
separately discuss the effect of the loss of spin order between the a6-planes at temperatures 
just above TV and spin misalignments inside these planes at higher temperatures. 

This work is organized as follows: in Section II we describe the system and the theoretical 
method, while in Section III we present the scattering formula in the Q symmetry, and its 
applications towards the description of the orbital and magnetic scattering in four special 
cases of spin arrangement. In Section IV, we present the corresponding orbital and magnetic 
spectra for Lao.sSri.sMnC^ and address the influence of the spin configuration, while in 
Section V we discuss our results in connection with the issue of the orbital spectrum of the 
paramagnetic phase. Section VI gives the conclusions of the work, while two appendices 
describe some calculational details. 



II. DESCRIPTION OF THE SYSTEM AND THEORETICAL APPROACH 

In the present work, we focus our attention on the half doped manganites with 
a low temperature CE antiferromagnetic structure. As an example, we consider the 
Lao.5Sri.5Mn04 compound, which displays a variety of phase o 1 ^ 1 ^ . The charge and or- 
bitally ordered phase sets on upon cooling at 230 K. The disproportionation of the Mn ions 
into the formal valences Mn 3+ and Mn 4+ is usually assumed, while recent exp er iment □ 19 1 20 
indicate that the charge difference in manganite materials is much smaller than 1. In the 
present atomic-based approach, we approximate the charge on the manganese ions by the 
integral values of 3+ and 4+. The four (three) 3d electrons of the Mn 3+ (the Mn 4+ ) are 
under the strong influence of the crystal field. In the presence of the cubic ligand field, 
the degenerate d orbitals split into triple degenerate t 2g and doubly degenerate e g levels. 
Cooperative Jahn- Teller distortions of the MnOg octahedra further reduce the symmetry of 
the crystal field to at the Mn 3+ site and thus lift the orbital degeneracy. t 2g level splits 
into d xy and doubly-degenerate d XZ)yz orbitals, and e g into orbitals of the d^ z 2_ r 2 and d x 2_ y 2 
type. Only one of the e g orbitals is occupied on the Mn 3+ ions, leading to orbital ordering, 
i.e. an alternate arrangement of orbitals oriented along the a and b axes on neighboring 
Mn 3+ A and B sites (see Fig. The antiferromagnetic phase, characterized as the CE type 
structure, exists under 115 K. In this phase, the spins are lying in the afr-plane 2 ^, displaying 
ferromagnetic ordering along spin-chains in the afe-plane (see Fig. H A~ — — A~ and 
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FIG. 1: The Mn02 plane of Lao.sSri^MnC^ representing charge, orbital and spin ordering for two 
different cases of spin orientation in the CE structure. Case I: the spin and crystal axis b enclose 
an angle a = 45° and Case II: a = 0°.The dashed lines denote the ferromagnetic spin chains. For 
reasons of clarity, the spins on the "Mn 4+ " ions are omitted. We assumed the x 2 — z 2 /y 2 — z 2 type 
of orbital ordering^, in the crystal coordinate system with x, y, and z along the a, b, and c axis, 
respectively. 



A + — B + — A + spin chains, "+" and "-" refer to the spin direction) and antiferromagnetic 
ordering between the spin chains both in the afr-plane and along the c-axis. 

Figure Q presents two possible cases of spin orientation in the CE structure, which we will 
consider in the present work. In Case I, the spins are parallel to the a + b axis, corresponding 
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FIG. 2: The Mn02 plane of Lao.sSri.sMnO^ representing charge, orbital and spin ordering for 
two hypothetical 90°-spin-canted structures. 



to an angle a = 45° between the spin direction and the crystal b axis. In Case II, instead, 
the spins are parallel to the b axis (a = 0°). 

Having in mind the paramagnetic phase and the effect of spin misorientations on the 
orbital spectrum, we will also consider two 90°-spin-canted structures, which are illustrated 
in Fig. |2J These structures, "A-type canted" and "B-type canted" , are built using local 
atomic configurations of Case II in Fig^ The A-type canted is obtained using the Mn 3+ 
A^ + '~' ) configurations (spin parallel to the orbital) and the B-type canted using the Mn 3+ 
B( +, ~) configurations (spin perpendicular to the orbital), with possible 90° rotations of these 
atomic building blocks. We note that, for the orbital scattering, only Mn 3+ sites need to 
be considered. Also, for the magnetic scattering of the CE structures, we will consider (as 
in the experiments^^) only superlattice reflections which selectively probe the Mn 3+ spins 
contribution. 

We based our calculations of the RXS spectra on atomic multiplet formulation in a 
crystal fiel d 22 ' 2 ? . Cowan's atomic multiplet progra m 24 ^ provides the Hartree-Fock values of 
the radial Coulomb (Slater) integrals F°' 2,4 (d,d), F 2 (p,d), G 1,3 (p,d) (direct and exchange 
contributions) and the spin-orbit interactions ((2p) and ((3d) for an isolated Mn 3+ ion at 
K. They are given in Table |l] It is known that the strengths of the Coulomb integrals in 
a solid are reduced due to the screening effects to about ~ 70-80% of the atomic values 2 ^. 
The largest of these integrals are nearly an order of magnitude greater than the crystal field 



splittings. This fact justifies our choice of the atomic multiplet approach in a crystal field. 



F 2 (d,d) F 4 (d,d) F 2 (p,d) G^d) G 3 (p,d) C(2p) C(3d) 
3d 4 11.415 7.148 - 0.046 

2p 5 3d 5 12.210 7.649 6.988 5.179 2.945 7.467 0.059 

TABLE I: The Hartree-Fock values for the ground and excited states of Mn 3+ given in eV. The 
p-shell spin-orbit parameter C(2p) has been increased by 9 % from the Hartree-Fock value to 
correspond to the experimental value^i. 

The treatment of the crystal field in this model is based on symmetry considerations. 
Starting from the transition matrix elements evaluated in the spherical point group, we lower 
the symmetry by the inclusion of a crystal field, following a "branching" of point groups. For 
example, the following "branching" can be generated: O3 D Oh D D D 2 h 3 Cffi For 
the calculation of the matrix elements in a given point group, the Wigner-Eckart theorem 
is used. In practice, this theorem is applied through a multiplication of the matrix element 
evaluated in the O3 point group by a chain of 2-j and 3-j coefficients for each branching. 
This is implemented in a crystal field code, which gives the energy spectrum and the dipole 
transition matrix elements necessary for the calculation of the scattering spectra. The values 
of the cubic (X 400 ) and tetragonal (X 220 ) crystal field parameters are obtained from a fit 
to the experimental RXS spectra. In the present case, we will use the values obtained in 
Ref. H for Lao. 5 Sr 1 . 5 Mn04. These values are X 400 = 5.6 eV and X 220 = 3.75 eV, and 
correspond to an x 2 — z 2 /y 2 — z 2 type of orbital ordering (with the x, y, and z axes along 
a, b and c, respectively). All the Slater integrals were scaled down to 75 % of their atomic 
value o 2 ^ 2 ? . Furthermore, in order to simulate the effects of the spin ordering on an isolated 
ion (superexchange and direct exchange interactions), we introduce a magnetic field acting on 
the spin of the atom. This field splits additionally the S = 2 quintet into S z = —2, — 1, . . . , 2 
levels. Inclusion of the magnetic field favors the S z = —2 level as the ground state of the 
atom. The value of the Zeeman energy is set to 0.02 eV. 



III. SCATTERING FORMULA 



In the first part of this section we describe atomic scattering for a general orientation of 
the spin axes with respect to the local crystal field axes. Thereafter, we apply the atomic 
scattering formula to find expressions for the orbital scattering of the four different spin 
structures considered and the magnetic scattering for the two different orientations of the 
spins in the antiferromagnetic CE-phase. 

A. Systematic expansion of the atomic scattering amplitude in invariants of de- 
creasing symmetry 

To calculate a diffraction spectrum, we first evaluate the atomic scattering from a single 
site and then add up the contributions of the atoms in the whole crystal. The type of order 
is selected by the scattering vector q which describes the periodicity of the orbital, magnetic 
or charge order. Depending on the relative orientation of the local crystal field axis and the 
spin direction, the local symmetry can change. In the case in which the crystal and spin 
axes are collinear, the site symmetry corresponds to the one of the crystal field. We use a 
tetragonal crystal field^, described by the D^h point group. If the spin direction, however, 
differs from the crystal field axis, the local symmetry can be as low as that of the Cj point 
group. 

Following the approach by Carra and Thole^i, we write the scattering amplitude as a 
linear combination of the product of pairs of tensors of increasing rank, which transform ac- 
cording to the irreducible representations of the spherical group SO3. To lower the symmetry 
to the Ci point group, we need to branch each irreducible representation into subgroup of 
irreducible representations, following the chain: 3 D Oh D -D4/1 3 C^h D C 2 h D C Only 
the totally symmetric representations contribute to the atomic scattering amplitude. If we 
stop the branching at the point group, we find for the atomic scattering amplitude the 
following expression: 

f^,c 4h = |A{(e} ■ e,)(F 1;1 + - i [(e* f x e,) ■ z] (F 1;1 - F_ ls _0+ 

(e* ■ z)(e; • z)(2F 0;0 - - 

where A represents a scattering coefficient^!, ej(e^) is the incoming (outcoming) photon 



polarization, and F m;m i is defined as: 



E 



<0|Jit| n )( n |ji # | > 



(2) 



E Q - E n + hw + iT/2 

where m and m' denote polarization states and are the electric dipole operators defined in 
spherical coordinates. |0) represents the ground state with energy Eq and \n) an intermediate 
state with energy E n . The photon energy is %u and V represents the broadening due to 
the core-hole lifetime. The local coordinate system is chosen in such a way that the z- 
axis corresponds to the highest local symmetry axis, which is also the quantization (spin) 
direction. Eq. ^ is known as the Hannon- Trammel formula^, which was derived for the 
spherical symmetry, and is applicable up to local symmetries as low as that of the point 
group. Lowering the symmetry to the C2/1 point group, the following expression for the 
atomic scattering amplitude is found: 

cE 

(3) 



/£,__, = &c 4h + |A{ [(e* f ■ y)(e 4 ■ y) - (e* • x)(e, • x)] (F_ 1;1 + F 1 ._ 1 ] 



(F_i.i-Fx._0}, 



(e}-x)(e l -y) + (e}-y)(e i -x) 

where x, y, and z are unit vectors corresponding to a local right Cartesian coordinate 
system, with the z-axis defined as above. Finally, assuming no symmetry operation on the 
site, except the inversion, we obtain four additional terms in the description of the atomic 
scattering amplitude in the Cj point group: 
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(Fl ; o + F_i ;0 + F 0; i + F 0; _i) + 

( — Fi ;0 — F 0; _i + F_i ; o + F 0; i) + 

(Fi;o + F_i ;0 — F 0; i — F 0; _i) + 

( — F 1;0 + F 0; _i + F_i;o — F 0; i)}, 

(4) 



where x, y, and z are unit vectors of any right Cartesian coordinate system, but it is 
convenient to choose, also in this case, the local z-axis along the spin direction. 



B. Orbital order scattering 

In tetragonal CE structures, the wave vectors probing the antiferro orbital ordering are: 
g°° = {^f 1 — f , + f , I) in units of (^, ^, ^), where n, m, and I are integers and a, c are 
the lattice constants. Taking into account the structure factors, the scattering amplitude for 
orbital wave vectors is proportional to the following combination of the atomic scattering 




FIG. 3: The scattering plane with the definition of the scattering and azimuthal angles 9 and </>. 
ki(kf) is the incident (outcoming) photon wavevector, q = kf — fcj. For the case of orbital order 
scattering with q°° =(\,\,0), the reference axes are: ui = ^=[110], u*2 = [001] and U3 = -7= [110] 

in the crystal coordinate system a, b, c. For the magnetic scattering, q AF ={\,—j,h) and the 
following reference axes are used: u*i = -4= [111], U2 = -4= [110] and 1I3 = -4= [112]. 



amplitudes at the A(+~) and B (+,_) sites: f A+ + f A ~ - f B+ - f B ~. In the soft x-ray 
experiments^"*^ the orbital order wave vector was taken as g°° =(|,|,0), and this is the 
wave vector we will consider for the polarization analysis presented in the next parts of this 
section. The equations will be presented using as reference system for the components of 
the atomic scattering tensors F m ^ m i the local coordinate system in which the z-axis remains 
in the spin direction, x lies in the a6-plane and the y-axis is perpendicular to the a6-plane. 

1. Case I: a = 45° 

For the case when a = 45°, it is useful to note that starting from the configuration at site 
A+, in Fig. Q and using symmetry operations (rotation by 180° about the c or (a + b) axis) 
one can obtain the configurations at the remaining three sites. This simplifies the numerical 
part of the calculation, because it suffices to perform an atomic calculation for just one site. 
The local symmetry of each site is described by the Q point group. From inspection of Eq. 
E]it can be seen that the only nonvanishing terms of the orbital order scattering are the ones 



including (e£ • x)(ej • z) and (e£ • z)(e, ■ x) polarization factors. Furthermore, the second 
term in the curly brackets in Eq. 0] vanishes, in our case, as the matrix elements in Eq. |21 
are real.— After evaluation of the polarization factor we obtain the following expression for 
the orbital order scattering amplitude: 

/£? = - v^2 cos 9 cos + <_+ + F% - F$+) , (5) 

in the case of a — > 7r' or 7r — > o 7 polarization. In Eq. El the scattering and azimuthal angles 
9 and are defined as shown in Fig. El using as reference axes ui = —^[110], u*2 = [001] and 
U3 = ^=[110] in the crystal coordinate system a, b, c. The scattering amplitude vanishes 
in the case of a — > a' or 7r — > n' polarization. The F term in parenthesis of Eq. El can be 
written in Cartesian coordinates as 2F XZ . 



2. Case II: a = 0° 



In the case when the angle a is different from 45°, the sites A and B are inequivalent, 
since the spin encloses two different angles, a a and with the local crystal field axis at 
sites A and B. Case II, with a = 0° in Fig. Q illustrates the situation with the highest 
local site symmetry, when the spin is collinear with the crystal field axis at site B (ag = 
0°), and perpendicular to the local crystal field axis at site A (a a = 90°). The local site 
symmetries for the B and A sites are C^ and C2/1 respectively. For the evaluation of orbital 
scattering, we add the contributions from the four sites in the same manner as we did in 
the previous case. But, this time, we need to do atomic calculations for two different sites, 
as there is no symmetry operation relating sites A and B. For the same reason, there will 
be no cancellation of contributions from various terms, which we had previously, and the 
expression describing orbital order scattering amplitude will have the following form for o 
polarized incident light: 
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where the upper expression in the curly brackets represents the a — > a' polarization part and 
the bottom one the a — > ir' contribution. The superscripts A and B denote the two different 
sites. The second term on the right hand side of Eq. |U] is proportional to F xx + F yy in the 
local Cartesian coordinates and the third one to F xx — F yy . The angles 9 and are defined 
in Fig 01 To simplify the polarization factors, we evaluate the polarization dependence for 
the azimuthal angle, using the value = 0°, as in the experiment (Ref. llOr ) . Taking this 
into account, Eq. |U]can be rewritten in a simpler form: 
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(7) 

with F^ B ^ + and F^- B ' + in the local Cartesian coordinates. We note the important fact 
that in the a = 0° case, the a — > a' polarization factor is nonvanishing, in contrast to the 
a = 45° case. 

In the case of n incident polarization, the scattering amplitude reads: 



f 
j i 



00,7Ti„ 

res 



Tjtcos 2 ^ - sin 2 9 sin 2 0) 
— | cos 0(sin 9 sin + cos i 



\ cos 2 9 — sin 2 6* (cos 2 (j> + \ sm 



\ cos 0(cos 9 + sin 9 sin (f>) 
\ cos 2 9 + sin 2 #(cos 2 - 
s cos 6(— 3 sin 9 sin + cos 



isin 2 




-^0;0 , 



+ 



F^U + Ftf - (F^U + F*i 



+ 



F% - + (F% + Ft, 



B+ 



(8) 

where the upper (lower) expression in the curly brackets corresponds to the ir — ► ir' (it — >• cr') 
polarization channel. Considering again = 0°, we see that all three terms on the right 
hand side of Eq. |H] contribute to the 7r — ► ir' scattering, while the expression for the it — > a' 
scattering is the same as in the case of the cx-polarized incident light. 



3. Canted structures 



The orbital scattering amplitude for the two canted structures in Fig. |2] can be described 
by the following expression with the F's in the local Cartesian coordinates: 



00,crinc 
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(9) 



assuming as above = 0°, and using the same convention as in Eqs. El and [7| for the 
polarization channels. For incident a polarization, the scattering takes place exclusively in 
the a — > ir' channel. 



C. Magnetic scattering 



We performed a similar analysis for the magnetic scattering of the antiferromagnetic CE 
structures, considering the wave vector q AF ={\,— which probes the antiferromagnetic 
order of the Mn 3+ atoms and was used in the experiment (Ref. llOl ). The corresponding 
magnetic scattering amplitude is proportional to: j af = f A+ + f B+ — f A ~ — f B ~ ■ In the 
case a = 45° the only contributing term is the second term on the right-hand side of Eq. d 
and, therefore, the antiferromagnetic scattering amplitude for incident cr-polarized light can 
be calculated from the expression: 

/4 F ' CTinc = 2 cos 6 cos (F A t ~ i^J-i) , (10) 



and in the case a = 0°: 



AF, trine 
res 



cos o[ cos H sin 0) H sin 6 

v 2 ^3 YJ 6 



- + - F^U) ■ (11) 



The angles 9 and have been defined as shown in Fig. 01 with the following definitions of 
the reference axes: u*i = -^[111], u 2 = ^[110] and u 3 = ^[112] in the crystal coordinate 
system a, b, c. There is no a — ► a' contribution and the polarization factor in Eqs. El and 
ITTlis evaluated for the o — >• ir' polarization. In the local Cartesian coordinates, the F— factor 
from Eqs. ITU1 and ITTI can be expressed as F xy — F yx . 



IV. INFLUENCE OF THE SPIN ORIENTATION ON THE SCATTERING SPEC- 
TRA 

In Fig. 01 we compare the calculated orbital order spectra of the CE structures with 
the two different spin orientations, a = 0° and a = 45°, displayed in Fig. ^ We present 
the spectra for a and 7r incident polarization evaluated using the equations derived in the 
previous section. In the case of the a incident polarization, panel a) displays the spectrum 
for the ct — > a' polarization, panel b) for the a — > ir' polarization, while the simulated exper- 
imental situation without a polarization analyzer for the outcoming photons is presented in 



panel c). Similarly, for the 7r incident polarization, we present the it — ■> tt' scattering channel 
on panel a), tt — > a' on panel b) and the average of the two in panel c). The polarization 
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FIG. 4: (Color online) Theoretical orbital order spectra of the CE structures with the two different 
spin orientations a = 45° (black line) and 0° (light green line) calculated for a and tt incident 
polarizations. Different panels represent the polarization channels: for the a incident a) a — ► a', 
b) a — ► 7r' and c) total, and for the tt incident a) tt — > tt' , b) tt — > a' and c) total. The experimental 
spectrum (dashed red line) for a incident polarization is reproduced from Ref. 



11 



dependence was calculated for the actual experimental angles: 9 = 62.9° and <fi = 0°—. We 
used T = 0.5 eV for the broadening due to the core-hole lifetime and the scattering inten- 
sity was convoluted with a Gaussian of width 0.1 eV to simulate the experimental (energy) 
resolution. 

The most significant difference between the spectra for the two angles a are the nonva- 
nishing a — > a' and 7r — > tt' contributions for the 0° alignment. For this angle and in the 
case of the a incident polarization, the scattering for the photon energies of the L3 edge 
is very small compared to the scattering at the L2 edge. For the a — » 71' polarization, the 
spectra for 0° and 45° alignments are identical, both in shape and magnitude. Similarly, for 
the 7r incident polarization, the ir — >• ir' contribution is much smaller than for the 7r — > a' 



scattering channel, the spectrum of which is again identical to the one obtained for a = 45°. 
This can be explained by considering a linear expansion of F A ( B ^ around the angle a = 45°. 
It can be shown (the details are given in Appendix B) that the first order corrections of four 
different F terms in Eq. [7| cancel exactly in the case of collinear spins, i. e. when the angle 
a a at the site A can be expressed as a a = 90° — leaving just the zeroth order term which 
equals the scattering amplitude in the case of 45°. In the situation where both polarization 
channels contribute with the same weight, [panel c)], for the a polarized incident light, the 
0° spectrum has smaller L3/L2 ratio because of the contribution from the a — > a' polariza- 
tion which is significant only at the L 2 edge. In the case of the 7r incident polarization, the 
branching ratio slightly increases for the unpolarized situation, as the n — > it' contribution at 
the L3 edge adds up to the main peak of the 7r — > a' polarization, while the L2 contribution 
causes appearance of a new small feature in the unpolarized spectrum. 

Considering only the situation c) (no outcoming polarization analysis), the orbital spec- 
tra calculated for the a = 0° and a = 45° structures show similar agreement with the 
experimental spectrum in Fig. HJj). In this situation, one cannot discriminate between 
the two cases. Very recently, however, experiments with a polarization analysis were per- 
formed on Lao.5Sri.5Mn04 (Ref. 0), and no detectable orbital scattering was observed for 
the a — > a' channel. Based on this observation, we conclude that the spin orientation in 
La .5Sr L5 MnO 4 should correspond to an angle a equal to, or very close to, the value of 45°, 
since a measurable a — * a 1 component^ should be present otherwise (see also Appendix A). 

The magnetic spectra for the two CE structures are compared in Fig. The main 
difference between the a = 0° and a = 45° cases is the decrease of the L 3 /L 2 branching ratio 
in the case a = 0°. Otherwise, the spectra have a similar shape. We note that a = 45° is also 
the value which was used in Ref. |ll| to examine the effect of the Jahn- Teller distortion on 
the orbital and magnetic spectra of Lao.sSri.sMnC^. The experimental orbital and magnetic 
scattering spectra and fits for La ^Sii 5 Mn0 4 have been presented and discussed in Ref. 



Our finding that a ~ 45° differs from the estimate of Staub et a/.—, who predicted 
a — 10 ± 5°, based on a fit focused on the contribution of minority ferromagnetic domains 
to the observed orbital scattering. We note, however, that the angle a in the minority 
ferromagnetic domains and in the majority CE domains may not be necessarily the same. 

To illustrate further the importance of taking into account the proper spin arrangement 




FIG. 5: (Color online) Theoretical magnetic order spectra of the CE structures for the two 
different angles enclosed by the spin and b axes: a =45° (black line) and 0° (light green line). The 
experimental spectrum (dashed red line) is reproduced from Ref. I 111 . 
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FIG. 6: (Color online) Theoretical spectra of the orbital scattering for the A-type (black line) and 
B-type (light green line) 90°-spin-canted structures, displayed in Fig. El 



(and all four different scattering sites), we calculated the orbital scattering from the canted 
spin structures, in Fig. El The orbital spectra for the two canted structures are displayed in 
Fig. IH1 It is clear that these two configurations yield spectra largely different from the exper- 
imental and theoretical spectra for a = 45°, and that the calculation of the orbital scattering 



in Lao.sSr! 5 Mn04, as shown in Fig. 01 should not be approximated by the consideration of 
only one site of high symmetry, A or B. At the same time, even from a quick inspection 
of Fig. IH1 one readily concludes that the two spectra are significantly different from each 
other. Yet, the only difference in their calculation was the spin arrangement. From this, 
we conclude that changing the spin arrangement can have a drastic influence on the orbital 
spectrum. This must be related, of course, in some way to the spin-orbit interaction. To bet- 
ter identify the origin of this effect, we performed calculations with either (i) the d spin-orbit 
parameter set to zero or (ii) the p-d interaction set to zero. Calculations (i) yielded only 
minute changes to the spectra. Calculations (ii), instead, yielded major changes, resulting 
in two virtually identical spectra for the A and B canted structures. Therefore, we conclude 
that the sensitivity to the spin direction comes from the combined effect of the spin-orbit 
interaction in the p-shell and the p-d Coulomb interaction, through which p-electrons feel 
the crystal field^. 

V. EFFECT OF DIFFERENT TYPES OF SPIN MISORIENTATIONS ON THE 
ORBITAL SPECTRUM IN THE PARAMAGNETIC PHASE 

With the increase of temperature T above T^v, spin ordering is first lost along the c- 
axis, while the in-plane correlations reduce gradually until the charge and orbital ordering 
temperature Tco= 230 K, at which they vanish completely^. Based on the results from 
the previous sections, we discuss the influence of various spin misorientations on the orbital 
spectrum in the paramagnetic phase for temperatures T^ < T <C Tco-, when there are 
still correlations inside the ah— plane. Randomness is expected to be mostly in the spin 
orientation between the planes, while the antiferromagnetic ordering inside the planes is 
assumed to be less influenced. In the present work we do not consider any change in 
the orbital configurations within the crystal. We calculate the scattering contributions at 
q =(\,\,0) from CE afr-planes with different spin orientations for the angles belonging to 
the class of a = 45° + n x 90° (n integer), which we defined as C(45°), and the class of 
a = 0° + n x 90°, C(0°), i.e. for a = 0°, 45°, 90°, 135°, 180° 225°, 270° and 315°. These 
planes are stacked in a random order, but assuming an equal proportion of each of them. 
In this case, the scattering amplitude is given by the average of the scattering amplitudes 
of the a = 45°, 135°, 0°, and 90° cases. The corresponding orbital spectrum is shown in 
panel a) of Fig. for the incident a polarization. We see that the spectrum does not change 
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FIG. 7: (Color online) Comparison of the theoretical orbital scattering spectra for a = 45° (black 
line), and for configurations with different types of spin misorientations in the paramagnetic case 
(light green line), described in the text. 

with respect to the spectrum obtained for the CE phase (a = 45°). This is because the 
a — > 7r' scattering amplitudes for a = 0°, 45°, 90°, and 135° are all equal, and with an equal 
amount of a = 0° and a = 90° layers, the atomic scattering amplitude for the o — > a' 
scattering component vanishes, reducing the spectrum to the a — > n' scattering component. 
Furthermore, assuming also a loss of coherence in the spin orientation between ferromagnetic 
chains within the a6-plane and a random occurrence of the angles belonging to the classes 
C(45°) and C(0°) for the ferromagnetic chains, the scattering amplitude can also be shown 
to be that corresponding to the average of the a = 0°, a = 90°, and a = 45° amplitudes. 

Hence, to produce a significant change in the orbital scattering, it is necessary to consider 
other changes (with a higher degree of misorientation) in the spin configurations, such as 
canted "ferromagnetic" chains. We find, however, that an equal amount of A and B-type 
canted configurations yields a spectrum identical to the one for the CE phase, and that 
it is necessary to consider a different amount of the two canted structures to change the 
spectrum. The orbital scattering spectrum for the case when 50 % of the sites corresponds 




to the configurations described with C(45°), 10 % to the A-type canted and 40 % to the 
B-type canted configurations is shown in b) panel of Fig. [7| The spectrum is characterized 
by the reduction of intensity of the first two peaks at the L 2 edge, as well as of the low- 
energy peak at the L 3 edge. The experimental data of Wilkins et a/.— indicate a large 
change of intensity at the L3 peak, a weaker one for the two low-energy peaks at the L2 
edge, and virtually no change for the high-energy L 2 peak. A qualitatively similar trend 
can be deduced from comparing the orbital spectra at 40 K and 160 K, measured by Staub 
et al&, for which the largest intensity reduction occurs at the L3 edge and the low-energy 
peaks at the L 2 edge, while the intensity of the highest-energy feature changes much less. A 
large change in the intensity of the main peak at the L 3 edge of the orbital spectrum at the 
Neel temperature could result from the magnetic scattering of the minority ferromagnetic 
domains, which would disappear above ~ Tn^. A similar effect, i.e. a significant reduction 
of the main peak at the L 3 edge and almost negligible changes of the other peaks, could also 
result from a decrease in the Jahn- Teller distortion^. However, the change in intensity of 
the other peaks is difficult to explain with the above effects, and Fig. [7]p indicates that they 
could be accounted for by a change in the spin configuration. 



VI. CONCLUSIONS 



In this work, we analyzed the dependence of the scattering spectra on the spin configura- 
tion in Lao.sSri.sMnC^for a fixed orbitally ordered configuration. We considered an arbitrary 
spin direction with respect to the local crystal field axis, yielding a low local symmetry, which 
necessitated a general scattering formula written for the C; point group. We generalized the 
Hannon- Trammel formula for the atomic scattering amplitude by a systematic expansion in 
invariants of decreasing symmetry up to Q and implemented it in the framework of atomic 
multiplet calculations in a crystal field. We discussed in detail the cases of four different 
spin structures. 

We first considered two antiferromagnetic CE structures with spins oriented either along 
the b axis or along an axis at an angle a = 45° relative to the b axis, within the afr-plane. 
The latter structure has the highest symmetry (all sites are related by symmetry operations) 
and has only the a — > 71' orbital scattering for a incident light which is confirmed by recent 
experiments which found no a — > a' scattering component. For the consideration of the 



CE structure with spins aligned with the lattice vector b (a = 0°), one needs to include an 
equal amount of sites which have spins either at 0° or at 90° with respect to the local crystal 
field axis in order to keep a collinear spin configuration. Orbital scattering in the a — > tt' 
channel for the cases a = 45° and a = 0° yields identical spectra, but in the latter case the 
scattering also takes place in the a — > tt' channel, in contrast to the experimental situation. 
Similarly, for the incident tt polarization, there is a nonvanishing tt — > tt' contribution, and 
the spectrum of the tt — > a' channel is identical to the one for a = 45°. For the magnetic 
scattering, the spectrum of the configuration with a = 45° displays a somewhat larger L2/L3 
ratio. From the comparison with experiment, we conclude that in the Lao.5Sr 15 Mn04 CE 
phase the angle between the spins and b axis is equal, or very close to, 45°. 

We then also examined two 90°-spin-canted structures obtained by considering only sites 
having the spin direction either at 0° or at 90° from the local crystal field axis. The spectra 
differ greatly in the two cases, and also differ from those obtained for the antiferromagnetic 
CE structures. From the consideration of these structures, we thus conclude that the spin 
configuration can have a drastic influence on the orbital spectrum. The largest differences 
caused by the spin configuration on the orbital spectra are observed in the presence of 
maximal disproportionation between spins parallel and perpendicular to the local crystal 
field axis. 

Finally, we also discussed the changes observed with temperature, in the orbital spectrum 
of the paramagnetic phase, considering the effect of various types of spin misorientations 
with respect to the CE structure. We first examined the effect of an incoherence in the spin 
orientation between the afr-planes and then also between the ferromagnetic chains within the 
afr-planes, considering ferromagnetic chains with spin directions corresponding to a = 0°, 
45°, 90°, 135°, 180°, 225°, 270°, and 315° in equal proportion. This was found to have 
a negligible effect on the orbital spectrum. We then examined the effect of a disorder in 
the spin direction within the ferromagnetic chains, considering the contributions of different 
types of 90°-spin-canted chains with more spins parallel than perpendicular to the local 
crystal field axis. This resulted in a spectrum with no change in the highest-energy peak, 
both at the L2 and L3 edges, but a reduced intensity of all other peaks. In view of explaining 
the measured temperature dependence of the orbital scattering intensity, the changes of the 
spectrum caused by the spin misorientations can complement the effects which would arise 
due to the presence of possible minority ferromagnetic domains or a reduction in the Jahn- 



Teller distortion, which would mostly affect the highest-energy peak at the L 3 edge. 



APPENDIX A: ORBITAL AND MAGNETIC SCATTERING FOR AN ARBI- 
TRARY ANGLE BETWEEN THE SPIN AND LATTICE AXES 

Previously, we discussed the scattering in the CE structure when the in-plane spin and 
lattice axes enclose angles of 0° or 45°. At this point, we present the expression for the 
orbital scattering amplitude in the case of an arbitrary angle a between the spin axis and 
lattice vector b for a-polarized incident light: 
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where we used /3 = 45° — a and the definitions of 9 and from Fig El We remind the reader 
that the F m ^ m i spherical tensor components are defined with respect to the spin axis, and 
that their values change with the angle a. We also note that there is an additional term in 
the scattering amplitude (the last term in the above equation), with respect to Eq. El which 
does not vanish in the general case. In Cartesian coordinates this term is proportional to 
F xz + F zx . Considering, e.g., the azimuthal angle = 0°, we note that there is a a — > a' 
scattering component for all angles < a < 90°, except for a = 45°. 

In the case of magnetic scattering, there are no additional terms (compared to Eqs. ITOl 
and lll|) contributing to the scattering. The scattering always takes place only in the a — > 71' 
channel. 



APPENDIX B: LINEAR EXPANSION OF F m m , AROUND a A[B) = 45° 

If we denote the angle between the spin and the local crystal field axis at site A (B) as 
ol a (olb), then the two angles in Fig. [TJ Case II, are related by: a a = 90° — as- We will 
expand the tensor F A m , defined in Eq. El as a function of the angle a a around 45° in the 
following way: 

*SU(«a) = <^(45°) + SF ^ aA) \ ( aA - 45°) + 0(Aa\). (Bl) 

The same expansion is valid for F^ m ,. Therefore, for the F A X and F zz from Eq. we can 
write: 

KMa) = 2^(45°) + 6 JM*A ^ (a A - 45°) + 0(Aa\), (B2) 
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where x is the coordinate along a and z the coordinate along b. We can also write: 
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using the fact that Aa A = -Aa B and that i^(45°) = F*(45°) and 5F ™ {aA) 
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, because the atomic scattering tensor F A (aA = 45° + Aa A ) at the A site is 

olb =45° 

equivalent to the atomic scattering tensor F b [olb = 45° + Aa A ) at the B site through a 
rotation of 180° about the a + b axis. In the intensity expression, Eq. F xx and are 
added up. From Eqs. IB2I and lB3l we see that the two first order terms will have the same 
magnitude and opposite signs. The same can be shown for the remaining terms from Eq. [?l 
F A and F xx , thus proving that the first order terms in the expression for the intensity all 
vanish, leaving only the zeroth order term equal to the 45° contribution. We also tested the 
assumption of a linear dependence of F m>m i on the angle around 45° by checking numerically 
the linear-order result: 

F*(a A = 0°) + F x A x (a A = 90°) = 2F A (a = 45°), (B4) 

and confirmed that the dependence F m>m t (a) is described well (within the numerical noise) 
with the expansion up to the first order. 
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When the p-d interaction is vanishing, and the d spin-orbit interaction is zero, the spin of the 
d electrons, Sd, and its rotations commute with the Hamiltonian. As a result, the scattering 
amplitude does not depend on the spin direction. On the other hand, when the p-d interaction is 
not vanishing, Sd and its rotations do not commute anymore with the Hamiltonian. Furthermore, 
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